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Aerodynamic  optimization  has  widely  become  a  issue  of  considerable  interest  to  determine  the  geometry 
of  an  aerodynamic  configuration  amidst  certain  design  constraints.  Aerodynamic  performance  is 
calculated  from  a  prescribed  geometric  shape,  which  is  often  performed  in  trial  and  error  method. 
Numerous  design  methods  are  available  for  the  aerodynamic  design  of  the  rotor.  The  goal  in  optimizing  is 
to  maximize  the  aerodynamic  efficiency  at  a  single  design  wind  speed.  However,  single-design  point 
methods  do  not  automatically  lead  to  the  optimum  design,  since  they  consider  only  one  point  in  the  total 
operational  range.  Moreover  they  do  not  implicitly  involve  considerations  on  loads  which  require  an 
experienced  designer.  The  aerodynamic  optimization  of  a  Horizontal  Axis  Wind  Turbine  is  a  complex 
method  characterized  by  numerous  trade-off  decisions  aimed  at  finding  the  optimum  overall  perfor¬ 
mance.  However  researcher  design  the  wind  turbine  is  an  enormous  ways  and  more  often  decision¬ 
making  is  very  difficult.  Commercial  turbines  have  been  derived  from  both  theoretical  and  empirical 
methods,  but  there  is  no  clear  evidence  on  which  of  these  is  optimal.  Turbine  blades  are  optimized  with 
the  aim  to  achieve  maximum  power  coefficient  for  the  given  blade  with  solidity,  ratio  of  coefficient  of 
drag  to  lift,  angle  of  attack  and  tip  speed  ratio.  In  this  article,  the  blade  element  theory  is  used  to  find  the 
optimum  value  analytically.  The  effect  of  power  coefficient  for  different  blade  angle,  tip  speed  ratio,  ratio 
of  coefficient  of  drag  and  coefficient  of  lift  and  blade  solidity  is  presented  and  the  optimized  set  value  is 
obtained. 
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1.  Introduction 

Wind  is  caused  due  to  uneven  heating  of  the  land  and  water  by 
the  sun  on  the  earth.  This  temperature  difference  induces  circula¬ 
tion  of  air  from  one  region  to  another.  Wind  energy  conversion 
system  is  one  with  low  investment  and  high  yield  power  genera¬ 
tion.  The  main  advantage  of  electricity  generation  from  wind 
energy  is  the  absence  of  harmful  emissions.  Thirst  on  energy,  the 
missile  price  increase  of  fossil  fuels,  hazardous  impact  on  environ¬ 
ment,  uncertain  supply  of  fuel  have  made  the  nation  to  rely  on 
renewable  energy  sources.  Variation  in  wind  velocity  round  the 
clock  must  be  extracted  with  the  available  wind  electric  converter 
and  hence  wind  electric  generators  must  be  capable  to  operate 
under  varying  speed  and  this  led  to  the  development  of  variable- 
speed  wind  turbines  nowadays.  The  power  produced  by  a  wind 
electric  generator  at  a  specific  site  depends  upon  the  mean  wind 
speed  at  the  site,  hub  height,  the  speed  characteristics  of  the  wind 
turbine  cut-in  speed,  rated  speed,  and  furling  wind  speeds. 

2.  Performance  and  design  aspects  of  horizontal  axis  wind 
turbine 


then  total  wind  power  available,  Pa, 

Pa  =  mw^m  ^yaffS  _  pAVm 

The  above  equation  shows  that  the  wind  power  varies  as  the 
cube  of  the  wind  velocity.  However  density  varies  with  pressure, 
temperature  and  relative  humidity.  Unfortunately,  the  total  wind 
energy  cannot  be  recovered  in  a  wind  turbine  because  the  output 
wind  velocity  cannot  be  reduced  to  zero,  otherwise  there  would  be 
no  flow  through  the  turbine  [2].  Let  Vt  be  the  inlet  wind  velocity 
and  V0  be  the  outlet  wind  velocity,  mw=p-  A  ■  Vave  be  the  mass  flow 
rate  with  average  velocity  (V,+V0)/2  and  the  power  recovered  from 
the  wind  ( Pout )  is  equal  to  the  rate  of  change  in  the  kinetic  energy. 

Pout  =  mw(Vf-V°)/2 
=  pAVave(Vj-V20y  2 
=  (pA)(Vi  +  V0)/2(V?-V20)/2 
=  (pA/4)(V?  +  vfV0-ViV2-V3o)/Vf 
Take  x  =  V„ /Vt 

Pout  =  Ci+X-X2-X3)Pa/2  (1) 


Designing  wind  turbines  to  achieve  satisfactory  levels  of 
performance  and  durability  should  have  deep  knowledge  in  the 
factors  affecting  wind  power,  aerodynamic  forces  acting  in  the 
turbine.  Energy  conservation,  pollution  prevention,  resource  effi¬ 
ciency,  systems  integration  and  life  cycle  costing  are  very  impor¬ 
tant  terms  for  sustainable  construction.  Designing  wind  turbine 
principles  includes:  (i)  minimizing  non-renewable  resource  con¬ 
sumption,  (ii)  enhancing  the  natural  environment  and  (iii)  elim¬ 
inating  or  minimizing  the  use  of  toxins,  thus  combining  energy 
efficiency  with  the  impact  of  materials  on  occupants  [1  ].  Therefore, 
possible  use  of  wind  energy  must  be  evaluated  in  terms  of  its 
impact  on  the  environment. 

2.1.  Power  available  in  the  wind 

The  three  factors  that  determine  the  output  from  a  wind 
energy  converter  are  wind  speed,  cross-section  of  wind  swept  by 
rotor,  and  overall  conversion  efficiency  of  the  rotor,  transmission 
system  and  generator.  Energy  available  in  wind  is  equal  to  the 
kinetic  energy  of  wind.  If  p  is  the  density  of  the  air  in  kg/m3,  A  is 
the  swept  area  in  m2,  and  Vm  is  the  mean  velocity  of  wind  in  m/s, 


Differentiating  (1)  with  respect  to  x  and  equating  to  zero,  we  get 
the  optimum  value  of  x  for  maximum  power  output. 

d(P0Ut)/dx  =  0=>1  -2x-3  x2  =  0 

Solving  the  quadratic  equation,  the  value  of  x=l/3. 
Substituting  the  value  of  x  in  (1),  we  get 

Pout  max.  —  0.593 Pa 


Thus  the  maximum  that  can  be  drawn  from  the  wind  system  is 
59.3%  of  the  total  wind  power  available,  which  is  called  Betz  limit 
in  aerodynamics. 

The  power  coefficient  Cp  is  defined  as 


Q>  = 


Ps 

(1/2  ).p.k.R2.U30 


(2) 


where,  Ps  is  the  shaft  power  output  in  Watts,  U0  is  the  upstream 
undisturbed  wind  speed  in  m/s.  The  power  performance  of  a  wind 
turbine  can  be  expressed  using  fixed  angular  velocity.  This  para¬ 
meter  is  defined  as 


(3) 
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Wind  turbines  have  various  Cp  values  depending  on  the  wind 
velocity.  Therefore,  their  efficiency  is  best  represented  by  a  Cp-k 
curve.  The  tip  speed  ratio,  X,  is  given  by 


where  X  is  tip  speed  ratio,  R  is  maximum  rotor  radius  (m),  a>  is 
rotor  speed  (rad/s)  and  V  is  wind  velocity  (m/s). 

The  available  wind  energy,  Ea  in  the  time  period  T  is  given 
by  [3] 

Ea  =  J^Padt  =  i pAjvldt  =  \pAVl  =  EasA  (5) 

where  Vm3  and  E ^  are  the  cubic  mean  wind  speed  and  the 
available  energy  flux  in  the  period  T. 

2.2.  Electrical  power  output 

The  power  in  the  wind  is  converted  into  mechanical  power 
with  power  coefficient  Cp,  with  generator  efficiency  >]g,  and 
mechanical  power  transmission  efficiency  t/m,  then  the  electrical 
power  output  Pe  is  given  as 

Pe  =  Cp.,lg.,lm.Pa  Watts  (6) 

Optimum  values  of  Cp,  rig  and  r/m  are  0.45,  0.9  and  0.95 
respectively  which  give  an  overall  efficiency  of  38%.  Actual  values 
will  probably  range  between  25  and  30%  [2]  which  may  vary  with 
wind  speed,  type  of  turbine  and  the  nature  of  load.  As  the  wind 
increases  from  a  low  value,  the  turbine  overcomes  all  mechanical 
and  electrical  losses  and  start  delivering  electrical  power  to  the 
load  at  cut  -  in  -  speed  V(  .  Rated  output  power  will  reach  at  rated 
wind  speed  VR,  above  which  constant  power  output  is  maintained. 
At  the  furling  speed  VF,  the  machine  is  shut  down  to  protect  it  from 
high  winds. 

The  efficiency  of  a  wind  turbine  is  usually  characterized  by  its 
power  coefficient  as  given  below.  Maximum  values  of  Cp  can  be 
0.5926  according  to  Betz  criteria  [4]. 


*~pe  —  2  3  v ' ) 

Omech-O  alternator -1^ /2)./>.R  X/ 

2.3.  Blade  design 

2.3.1.  Aerodynamic  design 

Blade  design  consists  of  aerodynamic  and  structural  design.  The 
challenge  is  designing  rotor  blades  for  different  applications  with 
optimised  weight  and  aerodynamic  performance.  The  major  appli¬ 
cations  of  rotor  blades  are:  (i).Stall  control  with  constant  speed,  (ii). 
Stall  control  with  variable  speed,  (iii).  Pitch  control  with  constant 
speed,  (iv).  Pitch  control  with  variable  speed.  The  usage  of  advanced 
design  tools,  production  technology  and  material  choice  are  the 
dominating  parameters.  As  the  power  curve  is  based  on  the  Cp  -  X 
characteristics  of  a  certain  rotor  blade,  the  main  parameters  for 
obtaining  an  optimum  power  curve  is  the  rotor  diameter,  rotational 
speed  and  pitch  angle  [5].  The  design  of  the  windmill  blade 
depends  on  the  following  parameters:  diameter  of  windmill  D, 
aerofoil  characteristics  and  the  number  of  blades,  Z  [6], 

2.3.2.  Production  technology 

With  composite  material,  production  technology  influences 
significantly  the  design.  In  the  field  of  rotor  blade  production, 
the  traditional  hand  lay-up  procedure  using  polyester  and/or 
epoxy  resin  as  matrix  material  together  with  glass  fibres  will  be 
substituted  in  the  near  future  by  more  advanced  technologies.  The 
application  of  preimpregnated  material  usually  suffer  from  too 
high  material  costs,  as  it  is  more  economic  to  use  the  raw 


materials  itself  in  the  production,  thus  using  resin  and  fibres.  An 
on-line  impregnation  technology  during  the  production  is  to  be 
used  for  an  adequate  production  of  unidirectional  stiffeners.  The 
usage  of  raw  materials,  such  as  filaments  of  glass  fibre,  is  generally 
said  cheaper  than  using  fabrics  with  different  lay-up  combina¬ 
tions.  Therefore,  a  compromise  has  to  be  found  between  the  usage 
of  the  different  raw  materials  and  its  consequences  in  design, 
especially  in  the  flexibility  of  the  structure. 

2.3.3.  Material  properties 

Basically,  there  exist  four  material  groups  used  for  rotor  blades: 
Epoxy  resin/glass  fibre,  Polyester  resin/glass  fibre,  Epoxy  resin/ 
wood,  Epoxy  resin/carbon-glass  fibres.  Further  improvements  in 
the  material  choice,  such  as  using  carbon  fibres  in  a  hybrid  system 
together  with  glass  fibres  has  been  used  for  rotor  blades  larger 
than  25  m,  as  a  sufficient  bending  stiffness  is  required  [7], 
Combined  with  an  optimised  structural  design  and  thick  profiles, 
it  is  also  possible  to  use  only  glass  fibres  for  rotor  blade  with  a 
length  of  30  m.  The  weight/strength  ratio  is  the  driving  parameter 
for  the  determination  of  the  optimum  material  and  lay-up  combi¬ 
nation.  Sandwich  structures  with  foams  are  necessary  for  the 
structural  stability.  On  the  other  hand,  material  damping  is  one  of 
the  mayor  issues  concerning  the  dynamic  behaviour  of  the 
complete  system  rotor  blades-wind  turbines,  especially  for 
epoxy/glass  fibre  and  polyester  resin/glass  fibre  systems.  The 
material  wood,  combined  with  epoxy  resin,  seems  theoretically 
to  have  excellent  performance.  Most  of  the  wind  turbine  blades 
are  made  of  fiberglass  reinforced  with  polyester  or  epoxy  resin. 
Small  turbine  blades  are  made  of  steel  or  aluminium,  but  the 
drawback  is  huge  weight.  Lighter  and  more  effective  blades 
decrease  material  requirements  for  other  wind  turbine  component 
making  overall  costs  to  be  lower.  Materials  with  lower  density 
such  as  fiber  aramid  (Technora)  have  higher  natural  frequencies 
and  bigger  deflection.  Blades  made  of  fiberglass  can  be  recon¬ 
structed  with  carbon  based  composites  to  reduce  mass  and 
increase  its  stiffness. 

Suggestions  for  increasing  performance  and  safety  of  windmill 
systems  listed  by  Onder  Ozgener  [4]  are  as  follows. 

•  Blades  can  be  made  of  epoxy-carbon  fiber  or  glass  fiber 
reinforced  plastics. 

•  To  produce  a  smooth  surface  a  steel  mold  can  be  used. 

•  A  long  and  narrow  airfoil  can  be  selected  having  larger  aspect 
ratio  than  the  classical  (short  and  wide  wing)  blade. 

•  Steel  blades  should  not  be  used  due  to  their  weight  and 
corrodibility. 

•  Lighting  protection  can  be  provided  for  GRP  epoxy-carbon  fiber 
blades. 

•  There  is  no  requirement  that  the  same  profile  should  be  used 
throughout  the  blade  length. 


2.3.4.  Weight 

The  blade  mass  is  one  of  the  most  important  parameters  for 
dynamic  loads  of  blade  and  wind  turbine.  The  aim  is  to  achieve  an 
optimum  between  a  low  weight  blade,  related  to  low-cost  produc¬ 
tion  and  a  high  performance.  The  blade  weight  can  be  reduced  by 
thick  profiles,  thus  increasing  the  moment  of  inertia  of  the  blade 
cross  section.  This  allows,  taking  into  account  the  material  elastic 
properties,  a  high  bending  stiffness. 

2.3.5.  Noise 

The  major  sources  of  noise  emission  of  rotor  blades  are  (i). 
Turbulence  in  flow  noise,  (ii).  Trailing  edge,  (iii).  Tip.  The  aero¬ 
dynamic  lay-out,  which  aims  an  optimum  aerodynamic 
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performance,  is  influenced  by  the  obligation  to  design  low-noise 
profiles  and  to  adapt  the  structural  lay-out,  especially  the  profile 
thickness  Furthermore,  dirt  on  the  blade  surface  contributes  to 
noise  emission. 

2.3.6.  Lightning  protection 

The  most  typical  damage  due  to  lightning  is  at  the  tip  region, 
where  the  increasing  temperature  leads  to  the  build-up  of  air 
pressure.  Therefore,  a  metal  alloy  receptor  is  integrated  in  the  tip 
region,  and  a  metal  stripe  inside  the  blades  transports  the  energy 
to  the  blade  root  connection.  From  there,  the  energy  is  transported 
to  earth  through  the  turbine  structure.  With  this  lightning  recep¬ 
tor,  it  is  only  necessary  to  repair  the  area  around  the  tip. 

2.4.  Wind  statistics 

Wind  is  a  highly  variable  power  source,  and  there  are  several 
methods  of  characterizing  this  variability.  The  most  common 
method  is  the  power  duration  curve  which  is  a  good  concept  but 
is  not  easily  used  to  select  Vc  and  Vr  for  a  given  wind  site,  which  is 
an  important  design  requirement.  Another  method  is  to  use  a 
statistical  representation,  particularly  a  Weibull  function  [2].  Local 
values  of  wind  velocity  should  be  3  m/s  or  higher,  and  the  wind 
should  be  steady,  to  produce  electricity  effectively. 

3.  Optimization  for  wind  turbine  rotors 

This  optimization  of  aeroturbine  focuses  on  the  development  of 
multi-disciplinary  optimization  algorithm  for  designing  of  hori¬ 
zontal  axis  wind  turbines  with  multiple  constraints.  The  aim  of  the 
optimization  process  is  to  optimum  potential  reduction  in  cost,  the 
optimum  specific  power  and  the  optimum  airfoil  characteristics. 
Design  variables  were  rotor  chord,  twist,  relative  thickness  and 
structural  shell  thickness  along  the  blades  with  the  tip  pitch  angle  [7], 

3.1.  Optimum  cost  of  energy 

To  compensate  the  reduction  in  annual  energy  production,  the 
swept  area  can  be  increased  to  gain  energy  yield,  without 
increasing  generator  size  and  design  fixed  loads  and  hence  total 
cost.  It  would  be  possible  to  constrain  the  energy  yield  to  a 
minimum  acceptable  value. 

3.2.  Optimum  rotor  geometry 

The  optimization  of  rotor  geometry  returns  smooth  shapes.  On 
reducing  the  chord,  the  blade  weight,  extreme  loads  and  fatigue 
loads  are  reduced  from  the  reduction  in  projected  blade  area.  The 
twist  in  the  root  region  is  of  minor  importance  to  the  power. 

3.3.  Optimum  airfoil  characteristics 

To  investigate  the  optimum  airfoil  characteristics,  the  lift  and 
drag  characteristics  should  be  considered  as  design  variables.  The 
high  at  the  root  is  often  studied  because  of  increased  production  at 
low  wind  speeds  before  rated  power.  A  reduction  in  chord  reduces 
both  blade  weight  and  extreme  loads,  but  should  be  counter¬ 
balanced  by  an  increase  in  the  CL  max  to  maintain  power. 

3.4.  Cost  of  energy  versus  specific  power 

Optimizations  were  done  with  different  constraints  on  the 
maximum  generator  power  to  investigate  the  variation  of  cost  of 
energy  with  the  specific  power.  Optimum  aerodynamic  efficiency 


at  some  design  wind  speed  is  closely  related  to  the  rotor  shape, 
however  the  efficiency  depends  on  the  number  of  rotor. 

3.5.  Blade  properties 

The  aerodynamic  profiles  of  wind  turbine  blades  have  crucial 
influence  on  aerodynamic  efficiency  of  wind  turbine.  When  blades 
of  length  more  than  45  m  are  used,  the  dynamic  behavior  of  the 
blade  should  be  taken  into  account,  and  the  position  and  shape  of 
the  spars  have  to  be  analyzed.  The  location  of  main  spar  together 
with  the  location  of  the  stiffness  ribs  will  have  the  biggest 
influence  on  the  bending  modes  of  the  blade.  The  twist  of  spars 
decides  about  pitch  of  principal  bending  axes. 

3.6.  Aerodynamic  loads 

Blade  Element  Momentum  (BEM)  method  is  used  for  the 
analysis  of  aerodynamic  loads.  It  is  an  iterative  method,  which 
assumes  the  value  of  axial  retardation  coefficient  ‘a’  to  be  zero  at 
the  beginning.  The  aerodynamic  loads  are  expressed  in  the 
following  formulas:  (c  is  the  chord  of  aerodynamic  profile) 


LiftL,  =  1. p.V2rel.c.CL  (8) 

Drag  ,D=^.p.V2rel.c.CD  (9) 

Thrust,  Fn  =  L.  cos  if  +  D.  sin  c/>  (10) 

Torque,  FT  =  L.  sin  </>-D.  cos  </>  (11) 

3.7.  Blade  material 


The  blade  is  made  of  composite  materials  with  more  than  one 
bonded  material  with  different  structural  properties  to  achieve  the 
combination  of  desirable  properties,  with  the  main  advantage  of 
high  ratio  of  stiffness  to  weight.  One  of  the  materials,  reinforcing 
phase  is  embedded  with  the  other  material,  matrix  phase.  The 
special  care  must  be  taken  in  defining  the  properties  and  orienta¬ 
tions  of  the  various  layers  since  each  layer  may  have  different 
orthotropic  material  properties.  Carbon  fiber  composites  allow  to 
less  blade  mass. 


4.  Composite  rotor  blades 

The  primary  objective  of  composite  rotor  blades  is  to  minimize 
the  blade  weight,  subject  to  frequency  and  auto  rotational  inertia 
constraints.  A  composite  is  a  structural  material  which  consists  of 
combining  two  or  more  constituents.  The  constituents  are  com¬ 
bined  at  microscopic  level  and  are  not  soluble  in  each  other.  One 
constituent  is  called  the  reinforcing  phase  and  the  one  which  is 
embedded  is  called  the  matrix.  The  reinforcing  phase  material  may 
be  in  forms  of  fibers,  particles  and  flakes.  The  matrix  phase 
materials  are  generally  continuous.  Strength  of  the  composite 
materials  depend  on  (i).Orientation  of  the  fiber  (ii).Type  and 
amount  of  the  fiber  present.  Fiber  orientation  in  each  layer  as 
well  as  stacking  sequence  of  the  plies  plays  a  major  role  in 
strength  and  modulus  of  the  composite  laminates.  Fibers  oriented 
in  one  direction  will  have  high  strength  and  stiffness  in  the 
direction  of  orientation. 

4.1.  Design  function  requirements 

The  design  functions  to  be  considered  are  as  follows 
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4.1.2.  Stiffness  and  strength 

A  combination  of  high  strength  and  stiffness  is  desirable 
because  of  the  vibration  from  the  natural  frequencies  in  the  air 
frame  and  the  periodic  loads  experienced  by  the  blade. 

4.1.2.  Weight 

The  most  important  fact  of  using  composite  material  is  con¬ 
siderable  weight  saving  which  is  determined  by  the  mass  moment 
of  inertia. 

4.1.3.  Safety 

Predictable  and  confidence  in  the  material  arises  only  with  the 
realistic  safety  margins,  to  maintain  safety  in  the  blades. 


4.1.4.  Impact  resistance 

The  blades  should  have  the  ability  to  resist  not  only  the  impact 
of  foreign  bodies  but  also  certain  level  of  mishandling  during 
servicing. 

4.1.5.  Erosion 

The  erosion  materials,  particles  in  the  air  such  as  dust,  sand  are 
very  abrasive  in  nature.  So  all  the  leading  edges  of  the  blades 
should  be  constructed  with  abrasive  resistant  materials. 

4.1.6.  Corrosion 

Corrosion  increases  the  safety  margins  and  decreases  the 
maintenance.  So  the  entire  part  of  the  blade  should  be  made  of 
corrosive  resistant  materials. 

4.1.7.  Cost 

The  main  design  optimization  of  composite  material  is  to 
satisfy  the  cost  requirement,  i.e.,  at  low  cost.  The  cost  includes 
low  initial  cost,  low  operating  cost  and  low  maintenance  cost. 

4.1.8.  Endurance 

Improving  the  survival  will  lead  to  high  reliability  and  less 
maintenance.  The  life  of  the  blades  has  important  implications  on 
operating  cost  and  must  be  maximized  to  ensure  economic 
viability. 

4.1.9.  De-icing 

A  facility  for  locally  heating  the  leading  edge  of  the  blade  is 
required  for  de-icing  purpose. 

4.1.10.  Lightning  strike  protection 

If  lightning  strikes  occur,  an  electrically  conductive  path  is 
required  along  the  blade  length  to  discharge  the  high  voltage. 


5.  Momentum  theory 

Wind  Turbine  extracts  kinetic  energy  from  the  wind.  Kinetic 
energy  in  the  wind  is  absorbed  by  wind  turbine  by  slowing  down 
the  wind.  If  it  is  assumed  that  mass  of  the  air  passing  through  the 
turbine  is  separated  from  the  mass  that  does  not  passing  through 
the  turbine,  then  the  separated  part  of  the  flow  field  remains  a 
long  stream  tube  which  lies  upstream  and  downstream  of  the 
turbine.  Fig.  1  shows  the  actuator  disk  model  of  a  wind  turbine. 

The  assumptions  made  in  this  theory  are  (i)  the  turbine  must 
be  a  horizontal  axis  configuration  such  that  an  average  stream 
tube  can  be  identified,  (ii)  The  portion  of  kinetic  energy  in  the 
swirl  component  of  velocity  in  the  wake  is  neglected  and  (iii)  the 
effect  of  the  radial  pressure  gradient  is  excluded.  The  upstream 
wind  velocity  V  is  decelerated  to  V(l-a)  at  the  turbine  disk  and  to 
V(1-2a)  in  the  wake  of  the  turbine.  Momentum  analysis  predicts 


the  axial  thrust  on  the  turbine  of  radius  R  to  be 

T  =  2;z-R2pV2a(l-a)(in  N)  (12) 

where  T  is  the  axial  thrust  on  the  wind  turbine  (N),  R  is  the  turbine 
radius  (in),  p  is  the  air  density  at  sea-level  at  standard  atmospheric 
conditions  (kg/m3),  Vis  the  wind  speed  (m/s),  ais  the  axial 
induction  factor. 

T 

Thrust  coefficient  CT  = - ..  =  4a(l-a)  (13) 

0.5  pnRfV2 

Mechanical  power  produced  by  the  turbine  is  given  as  P, 

P=27iR2pV3a(\-d)2{in\N)  (14) 

Wind  power  in  the  upstream  wind  covering  an  area  equal  to  rotor 
disk  Pw  (in  W), 

Pw  =  0.5pV3irR2  (15) 


Power  coefficient,  Cp  =  =  4a(l-a) 

Pw 


(16) 


Maximum  value  of  power  coefficient  is  at  a  =  1  /3,  substituting  the 
value  of  ‘a  in  Eq.  (16),  Cpmax  =  0.593 

Q 


Torque  coefficient,  Cq  = 


0.5  PnR2V2 


where  Q_  is  the  turbine  torque. 
CP 


Also,  C0  = 


(u/vy 


(17) 


6.  Blade  element  momentum  theory 

Blade-element  theory  helps  to  analyze  the  relationship 
between  the  individual  airfoil  properties  and  axial  induction 
factor,  power  produced  and  the  axial  thrust  of  the  turbine. 

The  elemental  torque  which  acts  on  all  blade  elements  in  an 
annular  ring  is 

dQ  =  0.5BcrpW2(CL  sin  </>-CDc os  </>)dr  (18) 

where  B  is  the  number  of  blades,  cis  the  chord  (m),  ris  the  radius  of 
blade  element  (m),  W  is  the  velocity  of  the  wind  relative  to  the 
airfoil  (m/s),  CL  is  the  lift  coefficient  and  CD  is  the  drag  coefficient, 
</>  is  the  flow  angle. 

The  sectional  lift  and  drag  coefficients  are  obtained  from 
empirical  airfoil  data  and  are  unique  functions  of  the  local  flow 
angle  of  attack  and  the  local  Reynolds  number  of  the  flow.  If  dland 
dD  are  the  lift  and  drag  forces  on  the  blade  element  respectively, 
then  lift  and  drag  coefficients  are  defined  as 

dL  =  CL(0.5pW2)c.dr  (19) 


dD  =  CD(0.5pW2)c.dr 


(20) 
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Power= torque  x  turbine  angular  velocity,  which  can  be  obtained 
by  integrating  Eq.  (18)  and  multiplying  the  same  with  angular 
velocity  of  turbine. 


rR 

P=0.5pBi2  /  cW2(Ci.sin  <p-CDc os  </>)dr 
Jo 

Similarly,  total  thrust  force  on  the  turbine  is  given  by 

rR 

T  =  0.5 pB  /  cW2(Ci.  cos  </>  +  CD  sin  <p)dr 
Jo 

From  the  Fig.  2, 
sin  $  =  V(l-a)/W 

cos  <p  =  (\  +a')r/W 

where  a'  is  the  tangential  induction  factor  =  co/2Q 
W  =  [V2(l-a)2  +  (1  +  a')2l22r2]05 

sin  a  =  sin  (4>-0) 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


For  flat  plate  airfoil,  CL  =  2 n  sin  a  ( 27 ) 

For  symmetric  airfoil  CL  =  2na  (28) 


For  circular  arc  airfoil,  CL  =  2 n[a  +  (2//c)]  (29) 

where /is  the  maximum  thickness  of  circular  arc  airfoil  (m) 

The  assumption  made  in  this  work  is  maximum  thickness  of  a 
circular-arc  airfoil  is  assumed  to  be  6%  of  the  chord,  hence  the  lift 
coefficient  equation  will  be 

CL  =  2n[a  +  (0.12)]  (30) 


Power  produced  by  the  wind  turbine  with  flate  symmetric 
airfoil  plate  can  be  obtained  by  substituting  Eq.  (23)-(27)  in  Eq. 
(10).  i.e„ 


P  =  jtpBQ^  J  V2.c.(l-a)2.c. cos  O.r.dr 
V.(l-a).(l  +  a').c.Q.  sin  0.r2.dr 
V.(l  — a).(l  +  a').c.£2(CD/Ci)..  cos  0.r2.dr 
+  Al  +  a')2.cJ22(CD/C,).sin  0.r3dr)  (31) 


Simplifying  assumptions  made  for  integrating  eq.  (31)  are 


i)  Uniform  distribution  of  upstream  wind  speed  along  the  blade 

ii)  Constant  chord  along  the  blade  (parallel  plan  form  blade) 

iii)  Constant  blade  angle  during  steady  state  operation 

iv)  Constant  drag-to-lift  coefficient  ratio 

v)  Uniform  distribution  of  axial  and  tangential  induction  factors 
along  the  blade. 


Introducing  the  definitions  of  solidity  and  power  coefficient 
and  integrating  the  above  Eq.  (20)  for  obtaining  power  coefficient, 

CP  =  w<7(l-a)(U/V){{  [l-a-(2/3)(l-a')(CD/CL)([//V)]cos  e 


(2/3)(l  +  a')(U/V)-0.5(CD/Ci.)(l  +  a')2(U/V)2 


(1-o) 


sin  0 


(32) 


Equate  Eqs.  (21)  and  (5)  to  get  the  axial  induction  factor, 
[nrr(U/V)(1—a) 


(2/3X1  +  a')(U/V) 


{  [l-a-(2/3)(l  +  cT)(Cd/Ci)(U/V)]  cos  0 
0.5(CD/Ct)(l  +  a')2(U/V)2 


(1-o) 


sin  6  (33) 


Similarly,  to  get  the  total  thrust  force  on  the  turbine,  substitute  the 
Eq.  (23) — (27)  in  Eq.  (22) 

T  =  npB  fvci-m+ayna*  drdr 

rR  rR 

-  /  (1  +  a')2cQ2  sin  0r2dr+  /  V2(1-a)2c(CD/CL)  cos  0dr 
Jo  Jo 

fR  1 

-/  V(l-a)(l  +  a')cX3(CD/CL) sin  Ordr  (34) 

Jo 

Integrating  Eq.  (34)  and  introducing  the  definition  of  solidity, 

T  =  pn2o(1-a)R2V2 ^  [0.5(1  +  a')(U/V)  +  (CD/Ct)(l-a)]  cos  6 


(1  +  a')2(U/V)2(l-a) 


0.5(CD/Ct)(l  +  a')(U/V ) 


sin  0 


(35) 


Equate  Eqs.  (35)  and  (12)  to  get  the  tangential  induction  factor. 


cos  0 


L  (U/V) 
-(1-a)(CD/CL)cos  0  + 

(l-a)2(U/V)2(l-a) 


0.5(U/V)cos  0 
2 


cos  0 


(U/V) 

0.5(1  +  a')(CD/Ci)(U/V) 


sin  0 


(36) 


Substituting  equations  for  circular  arc  airfoil 
P  =  rtpBi 2^1  V2 (1 -a)2 c cos  0rdr- J  V(1  — a)(l  +  a')ci2 sin  0r2dr 

rR 

-  /  V(l-a)(l  +  a')c£2(CD/CL) cos  0r2dr 
Jo 

+  J  (1  +  a')2cl22(CD/Ct)sin  flr3dr|  +  0.12ff/?6l2 
jy  V2(1-a)c(1 +Ar2)05rdr 
-Jo  V(1  — a)(l  +  a')cfl(CD/C[)(l  +/4r2)°-5r2dr| 

Where, 

A  =  (1  +a')2Q2/(1-a)2V2 


(37) 


CP  =  na(1  -a)(U/V)  { [1  -a-(2/3)(l  +  a')(CD/Cj(U/V)]  cos  0 

— [(2/3)(l  +  a’XU/V)- 0.5(1  +  a')2(CD/CL)(U/V)2/(1-a)sin  0]) 

+[0.2513ff(l-a)4/(U/V)(l  +  a')2]x(E15-l) 

-[0.1885(7(1  -a)3(CD/CL)E'5/(1  +  a')] 

+[0.0942(7(1 -a)3(CD/C1)£a5/(l  +  a')] 

+[0.0942a(l-a)4(CD/Ct)/(l  +  a')2(U/V )] 

xln|[(l  +  a')(U/V)/(l-a)]  +  £05|  (38) 
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where 

£=  1  +  [(1  +  a')2(t7/V)2/(l— a)2] 

a  =  [0.7854o-(U/V)/(l-a)]{[l-a-(2/3)(l  +  a')(CD/CL)(U /V)]  cos  9 
-[(2/3 K1  +  a')(U/V0-0.5(l  +  a')2(CD/Ct)(U/\/)2/(l-a)]sin  9} 
+[0.0628(7(1 -a)2/(l  +  a')2(U/V)](E15- 1) 

-[0.0471  ff(l  -a)(CD/Ci)/(l  +  a')]£15 

+[0.0236(7(1 -a)(CD/Ct)/(l  +  a')]E05 

+[0. 0236cr(l -a)2((CD/Ci.))/(l-a)2(L//V7)] 

xln|[(l  +  a')(fJ'/V)/(l  —a)]  +  £a5|  (39) 

T  =  /7/r2(7(l-a)R2V2{[0.5(l  +  a')(U/V )  +  (l-a)(CD/CL)] cos  0 
-[((l+a')2(U/V)2/3(l-a)) 

+0.5(1  +  a')(CD/CL){U/V)]  sin  0 } 

+[0.\257 npoRV3  (\-a)3 /(\  +  a')£, 2] 

+0. 1 885npoR2  V2(l  -a)2(CD/CL)E05 

+[0.1885;r/7(7W3(l-a)3(CD/CL)/(l  +  a'lflJE1  5 

-lxln|[(l  +  a')(t//V)/(l  —a)]  +  £05|  (40) 

a'  =  [1. 27324/(7(0// V)  cos  6»]{a-0.7854(7(U/V)  cos  6 
-1.5708<7(CD/Cj.Xl-a)cos  9 
—  1 ,5708cr[((l  +  a')2(U/V)2/3(l-a)) 

+0.5(CD/Ct)(l  +  a')(U/V)]  sin  0 
-[0.06285(7(1  -a)2  /(I  +  a’)(U/V)](Ev -5-l) 
-0.09425(7(l-a)(CD/Ct)£05 
-[0.09425(7(1 -a)2(CD/Ct)/(l-a)([//\7)] 

xln|[(l  +  a')(U/V)/(\-a)]  +  £°  5 1 }  (41) 

7.  Optimization  methodology 

The  theoretical  analysis  was  performed  to  investigate  the  effect 
and  dependence  of  the  various  parameters  in  the  wind  turbine 
rotor  geometry.  The  analysis  includes  the  recommended  values  at 
specified  operating  conditions  to  maximize  the  power  extracted 
by  the  wind  turbine  rotor.  In  this  work,  the  variation  in  the 
parameters  is  as  follows;  blade  angle  (9)  is  varied  between  0  and 
10°;  (CD/Ci)  is  varied  between  0  and  0.10;  rotor  solidity,  (a)  is 
varied  between  0.10  and  0.30;  tip  speed  ratio  (U/V)  is  varied 
between  2  and  14.  The  eqs.  (33)  and  (36),  since  they  are  coupled 


Fig.  4.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)  for  different  blade 
angle  G  with  (Cd/Q)  ratio  of  0.022  and  rotor  solidity  (<r)  as  0.03. 


Fig.  5.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade 
angle  6  with  (CD/Q)  ratio  of  0.022  and  rotor  solidity  (<r)  as  0.05. 


Fig.  3.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)  for  different  blade 
angle  0  with  (Cd/Cl)  ratio  of  0.022  and  rotor  solidity  (a)  as  0.01. 


Fig.  6.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade 
angle  Q  with  (CD/CL)  ratio  of  0.022  and  rotor  solidity  (<r)  as  0.07. 
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Fig.  7.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)for  different  blade  Fig.  10.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/l/)for  different  blade 
angle  B  with  (CD/CL)  ratio  of  0.022  and  rotor  solidity  (u)  as  0.1.  angle  6  with  (CD/CL)  ratio  of  0.022  and  rotor  solidity  (a)  as  0.25. 


Fig.  8.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (L//V) for  different  blade 
angle  6  with  (Cd/Cl)  ratio  of  0.022  and  rotor  solidity  (o-)  as  0.15. 


Fig.  11.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (Lf/VOfor  different  blade 
angle  0  with  (CD/CL)  ratio  of  0.023  and  rotor  solidity  (a)  as  0.01. 


Fig.  9.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 
angle  0  with  (Cq /Cl)  ratio  of  0.022  and  rotor  solidity  (o-)  as  0.2. 


Fig.  12.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade 
angle  6  with  (Cd/Cl)  ratio  of  0.023  and  rotor  solidity  (er)  as  0.03. 
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Fig.  13.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade 
angle  6  with  (Cd /Cl)  ratio  of  0.023  and  rotor  solidity  (<r)  as  0.05. 


Fig.  14.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 
angle  6  with  (CD/CL)  ratio  of  0.023  and  rotor  solidity  (a)  as  0.07. 


Fig.  16.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade 
angle  6  with  (CD/CL)  ratio  of  0.023  and  rotor  solidity  (<r)  as  0.15. 


Fig.  17.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade 
angle  0  with  (CD/Q)  ratio  of  0.023  and  rotor  solidity  (<r)  as  0.20. 


Fig.  15.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//V)for  different  blade  Fig.  18.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (l//\/)for  different  blade 
angle  6  with  (Cd/Cl)  ratio  of  0.023  and  rotor  solidity  (<r)  as  0.10.  angle  6  with  (Cd/Cl)  ratio  of  0.023  and  rotor  solidity  (<r)  as  0.25. 
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Fig.  19.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (L//V)for  different  blade  Fig.  22.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 

angle  B  with  ( CD/CL )  ratio  of  0.025  and  rotor  solidity  (u)  as  0.01.  angle  6  with  (CD/CL)  ratio  of  0.025  and  rotor  solidity  (a)  as  0.07. 


Fig.  20.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade  Fig.  23.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 

angle  B  with  (CD/CL)  ratio  of  0.025  and  rotor  solidity  (u)  as  0.03.  angle  8  with  (CD/CL)  ratio  of  0.025  and  rotor  solidity  (a)  as  0.10. 


Fig.  21.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade  Fig.  24.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 

angle  8  with  ( CD/CL )  ratio  of  0.025  and  rotor  solidity  (a)  as  0.05.  angle  8  with  (CD/CL)  ratio  of  0.025  and  rotor  solidity  (a)  as  0.15. 
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Fig.  25.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  [\3\V) for  different  blade  Fig.  2s.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (L//V) for  different  blade 
angle  e  with  (CD/Ct)  ratio  of  0.025  and  rotor  solidity  M  as  0.20.  angle  e  with  (Cd/Q)  ratio  of  0.026  and  rotor  solidity  (ff)  as  0.03. 


Fig.  26.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade  Fig.  29.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (L//V) for  different  blade 

angle  0  with  (Co /Cl)  ratio  of  0.025  and  rotor  solidity  (<r)  as  0.25  angle  9  with  (CD/CL)  ratio  of  0.026  and  rotor  solidity  (<x)  as  0.05 


Fig.  27.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade  Fig.  30.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade 
angle  9  with  (CD/CL)  ratio  of  0.026  and  rotor  solidity  (<r)  as  0.01.  angle  9  with  ( CD/CL )  ratio  of  0.026  and  rotor  solidity  (<r)  as  0.07. 
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Fig.  31.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade  Fig.  34.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 
angle  9  with  (CD/Ct)  ratio  of  0.026  and  rotor  solidity  (a)  as  0.10.  angle  9  with  (CD/Ct)  ratio  of  0.026  and  rotor  solidity  (a)  as  0.25. 


Fig.  32.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)for  different  blade  Fig.  35.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V) for  different  blade 
angle  6  with  ( CD/CL )  ratio  of  0.026  and  rotor  solidity  (<r)  as  0.15.  angle  9  with  ICD/CL)  ratio  of  0.028  and  rotor  solidity  (it)  as  0.01 


Fig.  33.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)for  different  blade  Fig.  36.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)for  different  blade 
angle  9  with  (CD/CL)  ratio  of  0.026  and  rotor  solidity  (a)  as  0.20.  angle  9  with  (CD/Ct)  ratio  of  0.028  and  rotor  solidity  (a)  as  0.03. 
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equations,  the  solution  of  those  equations  was  obtained  with  the 
help  of  Newton-Rapheson  two  variable  method  using  MATLAB. 
The  values  of  a,  a',(CD/Ct),  (U/V),  a  and  6  were  substituted  in  the 
equations  of  CP,  Ct,Cq  using  C  -  program.  The  optimized  or 
efficient  value  were  tabulated  by  finding  coefficient  of  variation. 

8.  Results 

Figs.  3-41  show  the  effect  of  power  coefficient  for  different 
blade  angle,  tip  speed  ratio,  ratio  of  coefficient  of  drag  and 
coefficient  of  lift  and  blade  solidity.  The  optimized  set  obtained 
is  at  blade  solidity  0.15,  angle  of  attack  5°,  tip  speed  ratio  8  and 
ratio  of  drag  coefficient  to  lift  coefficient  as  0.025. 

9.  Conclusion 

The  design  of  wind  energy  conversion  systems  is  a  very  complex 
task  and  requires  interdisciplinary  skills,  like  civil,  mechanical, 
electrical  and  electronics,  geography,  aerospace,  environmental 
etc.  An  attempt  has  been  made  to  discuss  the  important  design 
aspects  of  WECs.  The  prospering  future  in  wind  turbine  technology 
is  a  challenge  for  rotor  blade  design  in  order  to  enable  an  economic, 
reliable,  safe  and  maintenance  less  production  of  wind  energy.  This 


Fig.  37.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade 
angle  G  with  (Co /Cl)  ratio  of  0.028  and  rotor  solidity  (a)  as  0.07. 


paper  also  addresses  the  first  significant  step  by  defining  the  design 
parameters  of  developing  an  integrated  dynamic,  aerodynamic,  and 
structural  of  wind  turbine  blades  made  of  composite  materials. 


Fig.  39.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade 
angle  G  with  ( CD/CL )  ratio  of  0.028  and  rotor  solidity  (<x)  as  0.15. 


Fig.  40.  Effect  of  power  coefficient  Cp  with  tip  speed  ratio  (U/V) for  different  blade 
angle  G  with  (Cd/Cl)  ratio  of  0.028  and  rotor  solidity  (<r)  as  0.20. 


Fig.  38.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (L//V)for  different  blade  Fig.  41.  Effect  of  power  coefficient  CP  with  tip  speed  ratio  (U/V)for  different  blade 

angle  G  with  (Cd/Cl)  ratio  of  0.028  and  rotor  solidity  (<r)  as  0.10.  angle  G  with  (Cb/Ci)  ratio  of  0.028  and  rotor  solidity  (<r)  as  0.25. 
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Nomenclature 

Symbols 

A  swept  area  in  m2 

Pa  total  wind  power  available  in  Watts 

Pe  electrical  power  output  in  Watts 

mw  mass  flow  rate  of  the  wind  in  kg/s 

Vm  mean  velocity  of  wind  in  m/s 

Vi  inlet  wind  velocity  in  m/s 

V„  outlet  wind  velocity  in  m/s 

Vave  average  velocity  in  m/s 

Vc  cut-in-speed  in  m/s 

VR  rated  wind  speed  in  m/s 

VF  furling  speed  in  m/s 

Pout  power  recovered  from  the  wind  in  Watts 
Pout.max  maximum  power  that  can  be  drawn  from  the  wind 
in  Watts 

Ea  available  wind  energy 

Vm  mean  wind  speed 

Eos  available  energy  flux 

T  time  period 

Pe  electrical  power  output  in  Watts 

Cp  power  coefficient 

I  current  in  amps 

V  voltage  in  volts 

R  maximum  rotor  radius  in  m 

Ps  Shaft  power  output  in  Watts 

U0  upstream  undisturbed  wind  speed  in  m/s 
Ci  aerodynamic  lift  coefficient 

CD  aerodynamic  drag  coefficient 

CM  power  performance  of  a  wind  turbine 

c  chord  of  aerodynamic  profile 

L  lift  force 

D  drag  force 

Em  moment  force 

I  inclination  angle 

i  incidence  angle 

Greek  Symbols 
a  pitch  angle 


e 

angle  of  attack 

X 

tip  speed  ratio 

CO 

angular  rotor  speed  in  rad/s 

P 

Hellmann  coefficient 

>h 

generator  efficiency 

fjm 

mechanical  efficiency 

Va 

alternator  efficiency 

P 

density  of  air  in  kg/m3 

Abbreviations 

GRP 

Glass  fiber  Reinforced  Plastics 

NACA 

National  Advisory  Committee  of  Aeronautics 

HAWT 

Horizontal  Axis  Wind  Turbine 

VAWT 

Vertical  Axis  Wind  Turbine 

CSCF 

Constant  Speed  Constant  Frequency 

VSCF 

Variable  Speed  Constant  Frequency 

VSVF 

Variable  Speed  Variable  Frequency 

TSR 

Tip  Speed  Ratio 

0  &  M 

Operation  &  Maintenance  cost 

WECS 

Wind  Electric  Conversion  System 

IEC 

International  Electrotechnical  Commission 

SEIG 

Self-  Excited  Induction  Generators 
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